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Abstract 

Aims. We assess the importance of AGN outflows with respect to the metal enrichment of the intracluster medium (ICM) in galaxy clusters. 
Methods. We use combined ,/V-body and hydrodynamic simulations, along with a semi-numerical galaxy formation and evolution model. 
Using assumptions based on observations, we attribute outflows of metal-rich gas initiated by AGN activity to a certain fraction of our model 
galaxies. The gas is added to the model ICM, where the evolution of the metallicity distribution is calculated by the hydrodynamic simulations. 
For the parameters describing the AGN content of clusters and their outflow properties, we use the observationally most favorable values. 
Results. We find that AGNs have the potential to contribute significantly to the metal content of the ICM or even explain the complete 
abundance, which is typically ~0.5Zq in core regions. Furthermore, the metals end up being inhomogeneously distributed, in accordance with 
observations. 
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1. Introduction 

1.1. The ICM and Its Composition 

The main baryonic component of galaxy clusters is the in- 
tracluster medium (ICM), a hot (10 7 " 8 K), X-ray emitting 
plasma. Observati ons made with mod ern X-ray observatories 
like Chandra (e.g.lVikhlinin et all2005l) or XMM Newton (e.g. 
Tamu ra et al. 2004) lave revealed that the chemical composi- 
tion of the ICM is clearly non-primordial, with metal abun- 
dances of the order of ~Q.5Zq. Furthermore, metallicity maps 
are a vailable for the core regions of some very bright clusters 
(e.g. Sanders et all 120041 lHavakawa et all 120041 iDurret et al. 
2005). Apart from confirming the occurrence of high abun- 



Various groups h ave performed s imula tions on metal en 



dances, these also reveal that the metals are distributed inho- 
mogeneously. 

Since heavy elements can only be produced by means of 
stellar or explosive nucleosynthesis, it is obvious that they orig- 
inate mainly from cluster member galaxies. This demands for 
transfer processes which remove the metals from the galax- 
ies. Among the first suggested trans fer processes were ram- 
pressure snippi ng of interstellar gas dGunn & Gottlll972l) and 
galactic winds (IDe Young 1978). 
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richment in clusters. IDe Lucia et al.l (120041) used a combina 



tion of iV-body simulations and semi-analytic techniques to 
follow the formation, evolution and chemical enrichment of 
galaxies, including flows of gas and metals into and out of 
galaxies. They find that most of the metals currently in the 
ICM have been ejected at z > 1 and that massive galaxies 
are the most important conttibutors to the ICM metallicity. 



Nag ashima et al.l (|2005) took a similar approach with similar 



results. However, both g roups do not pr e dict th e distribution 
of the metals in a cluster. Tornatore et al. ( 20041) combined TV- 
body and SPH (smoothed particle hydrodynamics) simulations 
to create metallicity profiles of clusters, with a model that in- 
cludes star formation as well as supernovae to produce heavy 
elements, but did not distinguish between different metal trans- 
port processes. 

Our group distinguishes not only between different metal 
nansfer processes, but also between the ICM and gas that be- 
longs to galaxies. The enrichment efficiency of ram-pressure 
stripping, galac tic winds and merger-driven starburst s have 
been stu died bvlDomainko et al. ( 20061) . Kapferer et al. (2006) 
and ISchindler etaIT (l2005). Both processes were found to be 
able to significantly enrich the ICM, but they are not sufficient 
enough to explain the complete abundan ce. Galaxy-galaxy in 
teractions and intracluster supernovae (Ka pferer et alj 12005 
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Domainko et al . 2004) have also been found to be possible con- 
tributors to the metal content of the ICM. 

In this paper, we demonstrate that AGNs can also contribute 
significantly to the metallicity of the ICM. So far, investiga- 
tions on the impact of AGNs on galaxy clusters have been 
oriented almost exclusiv ely towards heating of the ICM (e.g. 
Siiacki & Springe! 2006), as this provides a possible solution 
to the cooling flow problem. In this work, however, we do not 
put the emphasis on energy transfer, but on metal transfer. 



1.2. AGNs and Outflows 

Observational evidence for wind-like outflows from AGNs ex- 
ists in the form of blueshifted absorption lines in UV and X- 
ray spectra. These lines are supposed to be created by matter 
moving away from the AGN, often with very high velocities 

of the order of thousands of km s _1 or even relativistic speeds 

I II — — k 

(e.g^ Chartas et al. 2003). Revi ews on t he subject can be found 

in fcrenshaw et all d2003h and IVeilleux et all d2005l) . The ab- 
sorption occurs at various distances from an AGN's central en- 
gine, from subkiloparsec to galactic scale, and in both radio- 
quiet and radio-loud AGNs. The inferred mass outflow rates 
are vague: ICrenshaw et al. I d2003l) concludes that they are com- 
parable to the mass accretion rates (which would likely leave 
a large range from gO.OOlMgyr -1 for low-luminosity AGNs 
to ~100M(r ) yr~ 1 fo r quasars and powerful radio galaxies) and 
Veille ux et al. (2005) propose the range 0.1 .. . IOMq yr 1 and 
suggest that ther e is a strong mass loa ding of the outflow by 
the galactic ISM. Morganti et al. (2005) assessed outflow rates 
ranging from ~1 to over 50Moyr _1 from the spectra of sev- 
eral radio-loud AGNs and estimated that the kinetic energy 
of the outflows might be sufficient to remove the gas from 
the bulges of the galaxies. A very distant (r ~ 28kpc), mas 



sive (~ 2 • 10 Mq) absorber has been found by lHamann et al 



(2001) in the vicinity of a radio-loud AGN, which indicates 
a mass outflow rate of ~10Moyr _1 provided that the mate- 
rial originates from the AGN and was expelled continuously 

np ositio n 
J d 19971) 

lower limit of ~Zq and says there is evidence for metallicities 
greater than IOZq. The possibility of super-solar ICM metal- 
li cities around AGNs h as been confirmed by the observations 
of llwasawa et alJd200lh . 



over an assumed lifetime of 10 yr. The compositio n of the 
outflows is likely to be rich in metals: lHaman nt dl997l) gives a 



The physical ori g in of these outflows is unclear 
dCrenshaw et al. 12003b iKrolikl 1 1999b : speculations include 
winds driven by hea t, radia tion pressure and magnetocentrifu- 
gal processes. Elvisl d2000b proposed an empirically derived, 
unifying geometrical model for the inner regions of AGNs 
that was assembled such as to explain all the absorption 
and emission features found in AGNs and their occurrence 
pattern. In this model, the absorption is ascribed to a warm 
(~10 6 K) and highly ionized medium, the "WHIM", which 
arises vertically from a narrow range of radii on the accretion 
disk in a funnel-shaped outflow. A self-consistent model which 
provides ready-to-use outflow properties does not exist yet, 
however. 



AGN jets in clusters are fou nd to be very exte nded, reach- 
ing -100 kpc into the ICM (e.g. lOwen et alll985b . Thus, they 
are an obvious candidate for the tr ansfer of metal-rich m atter 
into the ICM. The aforementioned iMorganti et alj ([2005), for 
example, ascribe their findings on interactions between radio 
jets and the surrounding material. The possibility that AGN 
entrainment is effic ient was confirm ed in early simulations of 
super-sonic jets bv lDe Young! d 1986b . who estimated that a to- 
tal of 1O 7 " 9 M0 of ambient material can be entrained by typical 
jets. This is interesting because an AGN must plough its way 
through metal-rich regions like the broad emiss ion line region 



The metallicity is probably very high there (Has inger et al 



l2002h . lBaldwin et alj (12003) found -15Z© for the BEL region 
in one particular quasar. The jets themselves, being outflows 
per se, are uninteresting with respect to metal enrichment, as 
they probably consi st either of pr otons and electrons or elec- 
trons and positrons dHirotanill2005b . 

Further evidence for galactic mass loss through AGNs is 
found in semi-analytic models impos ed on large-scale structure 
formation simulations: ICroton et alj d2006b find that AGNs are 
able to alter their host galaxies, thus inducing outflows of gas. 

2. Methods 

2.1. Numerical Methods 

We use combined A^-body and hydrodynamic techniques, to- 
gether with a semi-numerical galaxy formation code, to sim- 
ulate the different components of galaxy clusters. The dark 
matter A^-body simulations, which yield the gravitational po- 
tential of the dark m atter, are performed using a tree code 
dBarnes & Hull 1986b with constraint realizations of Gaussian 



rando m fields as initial c onditions ([ Hoffman & R ibak 1991 



van de Weygaert & Bertschingerl Il996[). The semi-numerica l 
galaxy formation and evolution code ( van Kampen et al] [l999) 
determines where galaxies form and provides galaxy proper- 
ties. Finally, the hydrodynamic code models the ICM. It uses 
the PPM ( piecewise parabolic metho d ) with a shock-capt uring 



scheme dColella & Woodward! 1 1984 iFryxell et all 1 1989) and 



computes the p roperties of the ICM in four nested, cubical grids 
(Ru fferlll 1992b . each of which is centered at the cluster center, 
with the largest having a volume of (20 Mpc) 3 and the smallest 
(2.5 Mpc) 3 . Each grid consists of 128 x 128 x 128 cells. Thus, 
we obtain the highest resolution in the cluster center. As initial 
condition for the gas, hydrostatic equilibrium is used. While 
the A^-body tree code and the semi-numerical galaxy forma- 
tion code start at z = 20, the hydrodynamic simulation starts 
at z = 1, thus covering only about 58% of the simulation time. 
We adopt a ACDM cosmology with Q A = 0.73, Q m = 0.27, 
<X8 — 0.93 and h - 0.1 . For more information on the methods, 
see lKapferer et al. I d2005b . 



2.2. AGN Outflows 

In the hydrodynamic simulations, outflows from AGNs are 
added to the model ICM at the respective position of the AGN 
host galaxy, altering the density, chemical composition, mo- 
mentum and energy of the ICM there. An outflow is either put 
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Table 1. Number of galaxies, average absolute value of the 
velocity and velocity dispersion (measured by the population 
standard deviation) of the three model clusters at z — 0. Only 
galaxies within r < 2 Mpc have been taken into account. 



cluster 


number of 


average velocity 


velocity dispersion 


galaxies 


[kms- 1 ] 


[kms- 1 ] 


1 


679 


1694 


828 


2 


396 


1078 


539 


3 


464 


1322 


623 



into the cell where the AGN host resides (one cell in the in- 
nermost grid has a volume of (19.5 kpc) 3 ) or it is distributed 
"jet-like" along a series of cells in the shape of bipolar out- 
flows, emanating from the launching AGN host up to 100 kpc 
to either side. 

To determine the AGN content and characterize the out- 
flows, a series of parameters has been introduced. In the ab- 
sence of a comprehensive theoretical understanding of AGN 
triggering and outflows, we have constrained these parame- 
ters using values from observations. The AGN fraction /agn 
is the fraction of galaxies that host AGNs with outflows . Multi- 
wavelength observations of clusters by iMartini et a l. (2006) 
have shown that the value for AGNs in general is about 5%. 
Assuming that 10% of all AGNs are radio-loud, we obtain 
/agn = 0.05 for outflows attributed to all AGNs (AGN winds) 
and /agn = 0.005 for outflows that are specific to radio-loud 
AGNs (entrainment by jets). We assume /agn to be fixed with 
time. Furthermore, we select AGN hosts rando mly among the 
cluste r members, relying upon the findings of iMiller & Owen 
(2002), who discovered that the radial distribution of AGNs is 
similar to that of cluster galaxies in general. 

The duty cycle or lifetim e tagn determines how long an 
outflow lasts. IMarti ni (2004) finds, collecting a broad variety 
of estimates, that the lifetime of AGNs is in the range 10 6 - 8 yr. 
He also notes that it is not clear whether the activity is episodic, 
with individual active periods much shorter than the total active 
lifetime. Whether outflows are present during the entire time is 
also arguable, of course. Nevertheless, we use the upper limit, 
tagn = 0.1 Gyr, as a standard value and refer to Sect. l4.2l for a 
discussion on the effect of shorter duty cycles. 

The mass outflow rate M out refers to the amount of mate- 
rial that an AGN galaxy loses to the ICM. As mentioned in 
Sect. 11.21 there are no good constraints on this quantity from 
observations, but outflows of the order of solar masses per year 
seem to be possible in at least some cases. The same holds true 
for the outfl ow metallicity Z nilt ; the value seems to be at least 
super-solar dlwasawa et alJEoOll) . We adopt 5Zq as a standard 
value. Note that the results for different metallicities are easily 
found, see Sect. l4.2l for details. 



Finally, there is the outflow temperature r out . We adopt the 
temperature of the WHIM (warm and highly ionize d medium) 
outflow (see Sect. II. 2t proposed by Elvis ( 2000h . which is 
10 6 K. 
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Figure 1. Number of galaxies within r < 2 Mpc with respect to 
the cluster center as a function of cosmic time. 
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Figure 2. Average absolute value of the velocity of all galaxies 
within r < 2 Mpc as a function of cosmic time. The humps 
indicate merging events. 



3. Properties of the Model Clusters 

Three different model clusters have been used in this work. 
Model cluster 1 is a rich cluster with an extraordinarily deep 
potential. Galaxies, detached or in small lumps, are attracted 
by its potential during the time interval in which the ICM is 
simulated. This and the formation of new galaxies raise the to- 
tal number of galaxies within r < 2 Mpc from 613 at z = 1 
to 679 at z = 0, see Fig.Q] This cluster also exhibits the high- 
est average velocities and the highest velocity dispersion, see 
Tab. [J 

Model cluster 2 differs considerably from model cluster 1 . 
It contains far less galaxies and is characterized by several 
merging events, which are visible as humps in the average ve- 
locity evolution curve shown in Fig. [2] At z = 1, the cluster 
consists mainly of two subclusters, separated by a distance of 
about 3 Mpc from each other. The subclusters approach each 
other and start to merge about 4.4 Gyr later. At z — 0, the 



4 



R. Moll et al.: Simulations of Metal Enrichment in Galaxy Clusters by AGN Outflows 



Table 2. Average galaxy number density n and average ICM 
density p within r < 1 Mpc with respect to the cluster center 
during the simulated time interval from z - 1 to z = (in 
runs without AGN outflows). The super- and subscripts give 
the greatest fluctuations from the mean. The ratio p/n, which is 
listed additionally, can be loosely interpreted as a "gas mass per 
galaxy". It characterizes the enrichment behavior of the clusters 
in our simulations: the ICM gets enriched more if p/n is low 
and vice versa. 



cluster 


n [Mpc 3 ] 


p [lO-^kgnr 3 ] 


£ [10 lu M o ] 


1 

2 
3 


96 6 +16 li 

45 r 26 - 2 

78 4+ 16 6 


22.2 +2 -* 
9.6!" 


4+0.4 

1 l +0 - 5 
i o+O.l 



Table 3. Parameter sets for the simulations presented in this 
section. From left to right: fraction of galaxies with AGN out- 
flows, duty cycle, outflow rate, outflow metallicity, outflow 
temperature and energy outflow rate (per AGN). The ejected 
material is distributed point-like in parameter set (a) and jet- 
like in parameter set (b), see text. 



/agn 


7" AGN 


M out 


Zout 


Tout 


^out 


[Gyr] 


[M yi->] 


[Z©] 


[K] 


[Js- 1 ] 


a 0.05 


0.1 


0.5 


5 


10 6 


6.0 ■ 10 32 


b 0.005 


0.1 


5 


5 


10" 


6.0 ■ 10 33 



galaxies settle down and form a single cluster with a total of 
396 galaxies within r < 2 Mpc. 

Model cluster 3 resembles model cluster 1 but contains less 
galaxies. With regard to its size and merging behavior, it can 
be considered as a standard galaxy cluster. The total number 
of galaxies within r < 2 Mpc rises from 449 at z = 1 to 464 at 
z — 0, albeit not mono tonic ally. This means that the present size 
of this cluster, as measured by the galaxy content, lies between 
cluster 1 and cluster 2. The same holds true for the average 
galaxy velocity and the velocity dispersion. 

The value p/n, by which we mean the average galaxy num- 
ber density divided by the average density of the ICM, charac- 
terizes the enrichment behavior of the clusters, as will be made 
clear in the following section. The values for our model clus- 
ters are listed in Tab. [2] In our model clusters, the depth of the 
potential determines the density of the ICM. The cluster with 
the deepest potential (1) is also richest in galaxies but it has 
nevertheless the largest p/n value. This of course relies on as- 
sumptions made in the simulations, such as the baryon fraction 
or criteria for when galaxies form. 

All our model clusters are "standard" clusters in the sense 
that they do not contain a cD galaxy or a cooling core. 

4. Results 

4.1. Metallicity Distributions 

The parameters for the simulations presented in this section are 
listed in Tab. [3] In the simulations (a) we assume that every 
AGN in the cluster has a moderate, wind-like outflow, whereas 
in the simulations (b) we assume that radio-loud AGNs ini- 



Table 4. Average metallicities from the 3d data (Z3&) and from 
the corresponding X-ray emission-weighted metallicity maps 
(Z map ) at z = 0. Zj,d is the average value of Z within a sphere of 
radius 1 Mpc, centered on the cluster center. Z map is the surface- 
brightness weighted average within a circle of radius 1 Mpc. 
Additionally, the mean number of AGNs within the mentioned 
region and the total mass of the metals that are directly released 
into the ICM there have been listed. 





Zw 


7 


mean # 


Zout-Mout,tot 




[Z©] 


[Z©] 


of AGNs 


[10 9 M o ] 


la 


0.027 


0.048 


22.6 


9.1 


lb 


0.030 


0.053 


2.3 


9.4 


2a 


0.082 


0.108 


10.0 


4.1 


2b 


0.097 


0.153 


1.1 


4.4 


3a 


0.048 


0.083 


18.0 


7.3 


3b 


0.050 


0.106 


1.9 


7.9 



tiate massive outflows through entrainment by extended jets. 
The material is initially distributed in just one cell in the for- 
mer case, and in several cells in the shape of bipolar jets in 
the latter, see Sect. 12.21 The specific values of the parameters 
are accounted for in Sect. 12.21 They have been tuned such that 
the total amount of metals released into the ICM is about equal 
in both cases, which facilitates a comparison. Forthwith, "2b" 
refers to model cluster 2 in simulation (b) etc. 

Figure [3] illustrates how the metals are added to the ICM, 
and how they are affected by the ICM dynamics in the case 
of the merging cluster 2 with jet-like outflows from radio-loud 
AGNs (b). At z = 1, the simulation starts with the ICM having 
Z = 0. Metal-rich material from AGNs is inserted at various 
positions. It has the same velocity as the AGN host galaxy and 
is accelerated along with the rest of the ICM towards the center 
of the cube, where the merging occurs. The jet-like shape in 
which the material is initially distributed is quickly disrupted 
by this movement. In the end, we find a blurred, but inhomo- 
geneous region of high metallicity in the cluster center. The 
simulation clearly shows that much of the metals that ended up 
in the center were ejected at much larger radii. 

The average metallicities in the cluster cores of the mod- 
els in all simulations can be found in Tab. H The emission 
weighted projection does not reveal the true metallicity; rather, 
it is approximately 60% higher. Although more than twice the 
amount of metals are ejected in the case of cluster 1 (the rich 
cluster) compared to cluster 2 (the poor, merging cluster), the 
latter has a metallicity which is on average more than three 
times as high. Cluster 3 lies in between. This correlates well 
with the values of p/n listed in Tab. [2] A cluster with a low 
p/n gets enriched more than a cluster with a high p/n and vice 
versa. 

The radial profiles of the 2D metallicity, depicted in Fig.|H 
show that most of the material ends up in the center of the clus- 
ters. Artificial X-ray emission weighted metallicity maps of the 
standard cluster 3 in both simulation scenarios are depicted in 
Figs. |5] and [6] The parameter set (b) was chosen such as to 
compensate a lower AGN fraction with a higher outflow rate. 
Nevertheless, the differences between (a) and (b) are quite dra- 
matic, as a comparison of Fig. [5] with Fig. [6] shows. The map 
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Figure 3. Distribution of the metals in simulation 2b. (a), (b), (c) and (f) show the evolution of the distribution in the innermost 
grid. The metals released into the ICM by the individual AGNs follow the merging process and get mixed quite thoroughly. 
Nevertheless, the final distribution is not homogeneous, (e) depicts slices through the final distribution, revealing that most 
metals reside in the core region, (d) shows what X-ray observations of the model cluster would ideally offer, if the model cluster 
was observed through an X-ray telescope with high resolution and high sensitivity. 
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0.5 1 1.5 2 

distance to X-ray maximum [Mpc] 



in Fig. [5] shows a pockmarked structure, whereas the map in 
Fig. |6]is dominated by single outflows that have not yet mixed 
with the surrounding material. 

The initial distribution of the outflowing material (point- 
like or jet-like) turned out to be only of minor consequence 
to the resulting metallicity distribution. The jet-like structures 
quickly (< 1 Gyr) dissolve. 

Considering the fact that we have only simulated about half 
of the time that might be at the disposal to AGNs for enriching 
the ICM, and keeping in mind that AGNs might have been more 
numerous before z - 1 , we see that both scenari os may explain 



the ob served metal abundances in the ICM (e.g. lSchmidt et al 
d2002l) find Z = 0.3 . . .0.6 in the core, r < 125 kpc, of the 



Perseus cluster). Also, in accordance with observations, the dis- 
tribution is invariably inhomogeneous. 



Figure 4. Radial metallicity profiles, created from the metallic- 
ity maps at z = 0. The bin width is 0.2 Mpc and the values in 
the bins have been weighted with the X-ray surface brightness. 
The profiles are slightly different for the two simulations (a) 
and (b). 



4.2. Parameter Variations 

Since the parameters describing the outflows and the outflow- 
capable AGN population of a cluster are not constrained very 
well, it is worthwile to make a parameter study. Of course, we 
know that the population of AGNs is very inhomogeneous. The 
way we used the parameters here is only a crude approximation 
to get a first notion of how they affect the resulting metallicity. 
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Figure 5. X-ray emission weighted metallicity map of model 
cluster 3 in simulation (a), where all AGNs are supposed to 
have wind-like outflows with O.SMqyi^ 1 , at z = 0. Fig. [6] 
shows the same for simulation (b). 




[Hpc] 



Figure 6. X-ray emission weighted metallicity map of model 
cluster 3 in simulation (b), where the radio-loud AGNs are sup- 
posed to eject metals into the ICM through entrainment by the 
jets, at z — 0. Although the total amount of metals is the same 
as in simulation (a), see Fig. [5] it shows a different structure 
with higher peak metallicities. 



Unsurprisingly, the average metallicity in a region covering 
the main cluster is linearly dependent on /agn, M oM and Z out . 
The difference between a low /agn and a high one is obvious 
from Figs. [5]and[6] 

Very massive outflows (>5Mq for standard runs, especially 
in clusters with a low pjn value and for the simulations (a), 
where the material was distributed into just one cell) can cre- 



ate features in the X-ray surface brightness and the temperature 
maps, where the ICM is both brighter and cooler than the sur- 
roundings. These features are caused by gas that is both cooler 
and denser than the surroundings, where material ejected by 
AGNs has not yet mixed thoroughly with the ambient medium. 
This is interesting because surface brightness maps from obser- 
vations are available at much higher resolution and for larger 
regions than metallicity maps. If lots of AGNs have massive 
outflows like that (e.g. 5% of the total galaxy population), both 
the surface brightness and temperature maps from the simula- 
tions reveal a flaky structure which is not compliant with ob- 
servations. Therefore, we discarded such models. 

Z ou t only affects the resulting metallicity values, not their 
distribution. This is to be understood as follows: a metallicity 
map resulting from a simulation where Z out = x can be repro- 
duced from a metallicity map from a similar simulation where 
Z ut = y by multiplying the latter with x/y. Given that r ou t is 
the same, a different Z ou t means a different E out (as the latter 
depends on the molecular weight), but the difference is minor. 
Thus, one can easily convert a map like that in Fig. |5]into one 
in which Z ou t is different by simply multiplying the values on 
the colorbar with the respective factor. 

Tagn determines how often the hosts of AGNs are changed. 
Thus, a shorter duty cycle means that the metal-rich gas from 
AGN outflows is dispersed more right from the start. This 
yields smoother metallicity maps. 

r ou t, together with M out and Z out , determines the energy out- 
flow rate E out . Putting a significant amount of energy into a cell 
of the model ICM results in vigorous convective movement, 
as a large portion of the energy is converted into kinetic en- 
ergy. The resulting metallicity map is much smoother and has a 
lower peak. For the simulations presented in the preceding sec- 
tions, this would happen if r ou t ~ 10 9 K. Below that, the value 
of the temperature is irrelevant to the final result. 

5. Conclusions and Outlook 

We find that, using a variety of assumptions based on obser- 
vations, outflows from AGNs can contribute significantly to 
the metal content of the ICM or even explain its complete 
abundance. Furthermore, the resulting metallicity distribution 
resulting from the simulations is always inhomogeneous, in 
agreement with observations. We showed that the allowable 
range for the AGN outflow rate can be constrained by com- 
paring the simulations with observations; general outflow rates 
> 5MQyr~' yield unrealistic results if the material is initially 
distributed in a volume of <(20kpc) 3 . 

However, AGNs are currently not understood well enough 
to definitely estimate their importance concerning metal en- 
richment. Apart from a better understanding of AGNs includ- 
ing triggering mechanisms and the origin of the apparent dif- 
ferences between them , the following questions need to be ad- 
dressed: What is the exact outflow rate of AGN outflows leav- 
ing the galaxy? How are outflows initially distributed in the 
ICM? Is the outflow continuous? How significant is the role 
of AGN jets? And how does all this depend on the properties 
of the host galaxy? A self consistent model of AGN outflows 
and jet simulations in a realistic, non-homogeneous environ- 
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ment would be helpful and will hopefully provide quantitative 
answers to some of these questions in the future. 

The cluster simulations are currently being improved by the 
implementation of new initial conditions for the hydrodynamic 
simulations that will allow us to start them at a higher redshift. 
Preliminary tests showed that this results into larger gas veloc- 
ities which leads to better intermixture between the material 
ejected by the galaxies and the ICM, especially in the outskirts 
of the model clusters. However, the inhomogeneity seems to be 
preserved. Future plans also include the discrimination of dif- 
ferent elements instead of only discriminating between hydro- 
gen&helium and metals. However, most important is certainly 
a better understanding of AGN outflows. 

Acknowledgements. This work was supported by the Austrian 
Science Foundation FWF (P15868, P18523 and P18416), 
Unilnfrastruktur 2005/2006, DFG (Zi 663/6.1), AUSTRIAN GRID, a 
University of Innsbruck scholarship and by ESO Mobilitatsstipendien 
of the Austrian Ministry of Science. We thank Chiara Ferrari for 
interesting discussions. 

References 

Baldwin, J. A., Hamann, F., Korista, K. T., et al. 2003, ApJ, 
583, 649 

Barnes, J. & Hut, P. 1986, Nature, 324, 446 
Chartas, G., Brandt, W. N., & Gallagher, S. C. 2003, ApJ, 595, 
85 

Colella, P. & Woodward, P. R. 1984, Journal of Computational 

Physics, 54, 174 
Crenshaw, D. M., Kraemer, S. B., & George, I. M. 2003, 

ARA&A,41, 117 
Croton, D. J., Springel, V., White, S. D. M., et al. 2006, 

MNRAS, 365, 11 
De Lucia, G., Kauffmann, G, & White, S. D. M. 2004, 

MNRAS, 349, 1101 
De Young, D. S. 1978, ApJ, 223, 47 
De Young, D. S. 1986, ApJ, 307, 62 

Domainko, W., Gitti, M., Schindler, S., & Kapferer, W. 2004, 
A&A, 425, L21 

Domainko, W., Mair, M., Kapferer, W., et al. 2006, A&A, 452, 
795 

Durret, F, Lima Neto, G. B., & Forman, W. 2005, A&A, 432, 
809 

Elvis, M. 2000, ApJ, 545, 63 

Fryxell, B., Miiller, E., & Arnett, D. 1989, MPA Report 449 
Gunn, J. E. & Gott, J. R. I. 1972, ApJ, 176, 1 
Hamann, F. 1997, ApJS, 109, 279 

Hamann, F. W., Barlow, T. A., Chaffee, F. C, Foltz, C. B., & 

Weymann, R. J. 2001, ApJ, 550, 142 
Hasinger, G, Schartel, N., & Komossa, S. 2002, ApJ, 573, L77 
Hayakawa, A., Furusho, T., Yamasaki, N. Y, Ishida, M., & 

Ohashi, T. 2004, PASJ, 56, 743 
Hirotani, K. 2005, ApJ, 619, 73 
Hoffman, Y. & Ribak, E. 1991, ApJ, 380, L5 
Iwasawa, K., Fabian, A. C, Allen, S. W., & Ettori, S. 2001, 

MNRAS, 328, L5 
Kapferer, W., Ferrari, C, Domainko, W., et al. 2006, A&A, 

447, 827 



Kapferer, W., Knapp, A., Schindler, S., Kimeswenger, S., & 
van Kampen, E. 2005, A&A, 438, 87 

Krolik, J. H. 1999, Active galactic nuclei : from the central 
black hole to the galactic environment (Active galactic nu- 
clei : from the central black hole to the galactic environment 
/Julian H. Krolik. Princeton, N. J. : Princeton University 
Press, cl999.) 

Martini, P. 2004, in Coevolution of Black Holes and Galaxies, 
ed. L. C. Ho (Cambridge: Cambridge Univ. Press), 170 

Martini, P., Kelson, D. D., Kim, E., Mulchaey, J. S., & Athey, 
A. A. 2006, ApJ, 644, 116 

Miller, N. A. & Owen, F. N. 2002, AJ, 124, 2453 

Morganti, R., Tadhunter, C. N., & Oosterloo, T. A. 2005, A&A, 
444, L9 

Nagashima, M., Lacey, C. G, Baugh, C. M., Frenk, C. S., & 

Cole, S. 2005, MNRAS, 358, 1247 
Owen, F. N., Odea, C. P., Inoue, M., & Eilek, J. A. 1985, ApJ, 

294, L85 
Ruffert, M. 1992, A&A, 265, 82 

Sanders, J. S., Fabian, A. C, Allen, S. W., & Schmidt, R. W. 

2004, MNRAS, 349, 952 
Schindler, S., Kapferer, W., Domainko, W., et al. 2005, A&A, 

435, L25 

Schmidt, R. W., Fabian, A. C, & Sanders, J. S. 2002, MNRAS, 
337,71 

Sijacki, D. & Springel, V. 2006, MNRAS, 366, 397 

Tamura, T., Kaastra, J. S., den Herder, J. W. A., Bleeker, 

J. A. M., & Peterson, J. R. 2004, A&A, 420, 135 
Tornatore, L., Borgani, S., Matteucci, F, Recchi, S., & Tozzi, 

P. 2004, MNRAS, 349, L19 
van de Weygaert, R. & Bertschinger, E. 1996, MNRAS, 281, 

84 

van Kampen, E., Jimenez, R., & Peacock, J. A. 1999, MNRAS, 
310, 43 

Veilleux, S., Cecil, G., & Bland-Hawthorn, J. 2005, ARA&A, 
43, 769 

Vikhlinin, A., Markevitch, M., Murray, S. S., et al. 2005, ApJ, 
628, 655 



